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ABSTRACT 

A two-level, linear, semiempirical theory for a failure criterion is described. The 
theory predicts the strength behavior of unidirectional filamentary composites under 
uniaxial and combined stress from basic constituent material properties and fabrication 
process considerations. Applications of the theory to several filament -nonmetallic ma- 
trix composites are presented and comparisons are made with experimental data. These 
results show good agreement between theory and experiment. Simple and combined 
strength envelopes are generated to illustrate the versatility of the theory and to point out 
problem areas in experimental work and design. 
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FAILURE CRITERIA FOR FILAMENTARY COMPOSITES- 
by Christos C. Chamis 
Lewis Research Center 

SUMMARY 

Failure criteria in the form of a two -level, linear, semiempirical theory to predict 
the strength behavior of unidirectional filamentary composites are presented. The 
first-level theory predicts the uniaxial strengths of the composite from its constituent 
material properties and fabrication process considerations. The second-level theory 
describes the strength behavior of the composite from its uniaxial (simple) strengths. 

The first -level theory is based on a modified rule of mixtures relation, on matrix-strain- 
magnification factors, and on maximum void effects. It considers both filaments and 
matrix as being generally orthotropic and reflects the particular fabrication process 
through the judicious incorporation of certain empirical factors. The second-level theory 
(combined-stress strength criterion) is based on a modified distortion energy principle. 

It is expressed as an interaction equation with a coefficient depending on the composite’s 
elastic properties. This coefficient is to be modified for theory -experiment correlation. 
The criterion is applicable to materials exhibiting different magnitudes in tensile and 
compressive strengths as well as to isotropic materials. 

The physical bases and the mathematical formalisms leading to the two-level theory 
are described. Suggested experimental techniques to measure the simple strengths and 
to evaluate the correlation coefficients (empirical factors) are illustrated. The need for 
controlled experimental data and complete test records is emphasized. Application of 
the two-level theory to several multilayered filament -matrix composites is presented to 
illustrate its application and to compare it with experimental data. These results indi- 
cate that the two -level theory predicts the composite strength behavior reasonably well. 
Composite simple and combined strength envelopes are generated to Illustrate the ver- 
satility of the theory and to point out problem areas in experimental work and design. 


* A portion of this work was supported by the Advanced Research Projects Agency, 
Department of Defense, through Grant no. AF 33(615) -3110, administrated by the Air 
Force Materials Laboratory while the author was a member of the Engineering Design 
Center, Case Western Reserve University, Cleveland, Ohio. 



INTRODUCTION 


Experimental observations have shown that unidirectional fiber composites exhibit 
five primary failure stresses, also referred to as limit stresses and simple strengths. 
These stresses result when the composites are subjected to failure under uniaxial loading 
in their plane. The failure stresses are identified individually as follows: (1) longitudi- 
nal tensile (2) longitudinal compressive (3) transverse tensile (S^22T^’ 

(4) transverse compressive (S^22C^’ intralaminar (inplane) shear (S^j^2s^‘ 

stresses are illustrated in figure 1 . (See also ref. 1, ch. 2.) Experimental observations 



®U2 1 hl2S, 

^I22ci°m^h22 cs- 49307 

Figure 1. - Unidirectional filamentary composite (geometry and simple 
strength definitions). 


also show that unidirectional fiber composites under combined loading fail at stress con- 
siderably different from those under simple loading (ref. 1, ch. 3). Therefore, five 
simple strengths and a combined-stress strength criterion are needed to describe the 
strength (limit) behavior of a unidirectional filamentary composite (UFC). Three alter- 
natives are possible in determining the simple strengths and the combined-stress strength 
criterion: (1) empirical, (2) theoretical, and (3) semiempirical. 

The empirical approach evaluates the simple strengths and the combined -stress 
strength criterion by measurement of specific filament-matrix systems. This approach, 
though effective if only one system is considered, is economically prohibitive when fea- 
sibility and trade-off studies are needed in considering several filament-matrix systems 
with a range of filament volume content for each system . The case for combined loading 
is far more complex. This requires complex experimental setups and many experiments 
to establish a reasonable failure criterion for a given load envelope. 

The theoretical approach is based on some mathematical model of the physical 
makeup of a UFC . Examination of its physical makeup reveals that its strength depends 
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on the properties of its constituents and the particular fabrication process. What is 
needed, then, is a mathematical formalism to relate the strength of the UFC to the prop- 
erties of its constituents and to the particular fabrication process. This mathematical 
formalism can be constructed with the aid of a two -level theory. The first level is a 
theory to predict the simple strengths of the UFC from constituent properties and fabri- 
cation process considerations. The second level is a theory to predict the onset of fail- 
ure (limiting condition) of the UFC from either the simple strengths (predicted in the 
first level) or from certain measured values. The failure criteria at both levels, to be 
effective and useful for design practice, must account directly or indirectly for the fol- 
lowing desirable features: 

(1) They must have some theoretical basis. 

(2) They should reflect the particular fabrication process (void content, size, and 
distribution; filament spacing nonuniformities and misalinement; differences in bulk 
and in-situ constituent properties; imperfect interface bond; residual stress, etc.). 

(3) They must be applicable to both isotropic and anisotropic filaments and matrices . 

(4) The resulting equations should be relatively simple. 

(5) The theories must be experimentally substantiated at both levels. 

Various theories have been proposed in the literature. Representative first-level 
theories Include: rule of mixtures (filament strength limited), netting analysis (ref. 2), 
and statistical (refs. 3 and 4) for longitudinal tensile strength mixtures 

(matrix strength limited) (ref. 5), filament microbuckling (ref. 6), panel buckling 
(ref. 7), and constituent debonding with intralaminar shear (ref. 8) for longitudinal com- 
pressive strength (S^jjq); and matrix stress concentration factors (refs. 9 and 10) and ma- 
trix strain magnification factors (refs. 11 and 12) for transverse tensile (S^22T^’ trans- 
verse compressive (S^22C^’ intralaminar shear (S^j2S^ strengths. Representative 
second-level theories include: statistical (ref. 13), curve -fit quadratic (refs. 14 to 16), 
special ellipsoids (refs. 16 and 17), distortion energy (refs. 18 and 19), maximum strain 
(ref. 20, combinations (ref. 21), and fracture mechanics (ref . 22). Many of these the- 
ories are extensively discussed in chapters 2 and 3 of reference 1 and chapter 3 of ref- 
erence 16. Here, suffice it to say that all these theories are deficient in either one or 
more of the desirable features delineated previously. 

The semiempirical approach has not been explored to its fullest extent. In this re- 
port a two-level, semiempirical theory is described which meets (directly or indirectly) 
all the desirable features delineated previously. The theory at both levels is linear and 
developed primarily for nonmetallic composites. The basic hypothesis of this approach 
is that variables which cannot be accounted for directly are incorporated indirectly 
through the judicious .introduction of theory -experiment correlation factors. The physical 
bases and justifications for these factors are discussed in this report. 
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SYMBOLS 


Ililllllli^^HWII I III ^ 

I 

] 

A array arrangement, distribution i 

;! 

A parameter, eq. (A9) 

a^,a£ constants, eq. (4) 

d diameter 

E loi^tudinal modulus ® 

F failure function 

G shear modulus 

K,K' elastic and correlation coefficients, respectively, eq. (12) 
k,k apparent and actual volume ratio, respectively 

N number of filaments or voids 

Njj, Ny applied loads 
p parameter, eq. (9) 

S simple strength, failure or limit stress 

t thickness 

/3 theory-ejq)eriment correlation factor 

e strain 

0 angle between load and filament directions 

V Poisson’s ratio 

a stress 

(p strain-magnification factor 

M 

Subscripts: 

B interface bonding 

C compression 

D debonding 

f filament property 

1 ply property 

m matrix property 

p limitii^ property 
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R 

S 

T 

V 


x,y,z 

1,2,3 

Q!,i3 


residual stress 
shear 
tension 
void 

load axes 

material axes (the 1-axis coincides with the filament direction) 
T or C tension or compression 


UNIAXIAL SIMPLE STRENGTHS 


The physical makeup of the UFC suggests. that its simple strength will be related 
to its constituent properties and to the fabrication processes as follows: 


= f (k,d,N,A)f 


(E,i',G,S, e ) , 

Pf,m 


S. 


B’ '^R 


( 1 ) 


where in equation (1) denotes the UFC simple strength; (k, d, N, A)j ^ denotes 

volume content, size, number, and distribution of filaments and voids; k^^ denotes the 

volume content of the matrix; (E, t-, G, S, represents the elastic and strength 

^ f, m 

properties of the filaments and matrix; and Sg and Sg denote interface bond strength 
and residual stress, respectively. The void content, the bond strength, and the residual 
stresses are dependent on the filament surface treatment. They also depend on various 
matrix additives, hardeners, temperature, and pressure during fabrication and on the 
fabrication method of making the UFC. As can be seen, the list of variables on which 
the UFC simple strei^ths depend is quite loi^. 

Evaluation of the function in equation (1) presents a formidable task which requires 
sophisticated statistical methods and a large number of experiments. This results in 
complex mathematical expressions not readily amenable to use in design. Variables 
such as void size and distribution, filament spacing nonuniformity, interface bond 
strength, and residual stresses are influenced by the particular fabrication process. 

If it is assumed that the particular fabrication process remains approximately invariant, 
then it is reasonable to group all these variables into theory -experiment correlation 
factors. The concept just stated simplifies the derivations and resulting expressions of 
the simple strengths and yet retains all the essential parts of equation (1) . The details 
are described in the following sections. 
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Longitudinal Tensile Failure Stress 


The expression describing the longitudinal tensile strength is a modified rule of 
mixtures equation of the form 


®Z11T 


- Sj,p 





( 2 ) 


where is the filament bundle strength (or single filament strength for monofilament 
composites); /Sjrp and are the theory -experiment correlation factors, which ac- 

count for the particular fabrication process; and are actual filament and matrix 
volume contents and are defined in the appendix; and and E^jj represent the in- 

situ longitudinal moduli of the matrix and filament, respectively. 

Several important points should be noted: (1) equation (2) is linear in kf if % 
and are independent of kj. The filament bundle strength (single filament for mono- 
filament composites) controls (2) When ®Z11T insensi- 

tive to iSjjjrp- (This is the basic hypothesis in the netting analysis.) (3) The coefficients 
/Sjrj. and should be relatively independent of the filament volume content for a fixed 

fabrication system. (4) The proximity to unity of these coefficients is a measure of the 
validity of the rule of mixtures and of the relative insensitivity of to the fabrication 

process. 


Longitudinal Compressive Failure Stress 


Using the rule of mixtures, the longitudinal compressive strength is related to con- 
stituent properties by the following equation: 


®mc 


= S 


mC 


^mC^m %C^f 



(3) 


Under compressive loading, it is also possible that a UFC will fail by a combination 
of debonding and intralaminar shear (ref. 8). This condition is approximated by the 
expression 


®aiD ~ ^1®Z12S ^2 


(4) 
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where is the matrix compressive strength, and are theory -experiment 

correlation factors analogous to those for a^ and ag are empirical curve-fit 

parameters, and S^j^2s UFC intralaminar shear strength, which will be defined 

subsequently. The remaining variables have already been defined. Equation (3) is the 
modified rule of mixtures matrix-strength limited but equation (4) limits by a 

combination of constituent debonding and intralaminar shear strength. The original 
forms of these equations were proposed by Fried (refs. 5 and 8). The equations were 
modified in reference 1 to incorporate several of the desirable features delineated in the 
INTRODUCTION. Equation (4) evolved from the experimental work of reference 8 where 
it was discovered that the longitudinal compressive failure stress depends on the intra- 
laminar shear strength. The curve -fit parameters a^ and ag in equation (4) are eval- 
uated as is described in references 1 and 8. It is possible that they would remain the 
same for various filament but only one -matrix systems. Though it is suspected that this 
might be the case, these coefficients should be evaluated for each filament-matrix sys- 
tem for reliable predictions. The important point to be noted in equation (3) is that 
S^llC is very sensitive to since important point to be noted 

in equation (4) is that the various fabrication process effects are introduced through 
S^ 12 s> which is defined subsequently. It is suggested as a conservative measure that 
S^llC be taken as the smaller of the two values computed from equations (3) and (4) or 
in equation form 


(5) 



r 

r 

- 


S^llC = MIN^ 

SmC 

^mC^m %C^f 1 

^%ii\ 

fmnj_ 

, (aiS^j2S ^2^^ 


Transverse Tensile Failure Stress Sj22j 

The governing equation for the transverse tensile strength is based on the hypothesis 
that S^22T limited by the allowable tensile strain in the matrix (ref. 1, section 2.5). 
In equation form this condition is expressed by 


h22T ~ ^^^mpT^ ^122 


( 6 ) 


where is the allowable matrix tensile strain defined in figure 2(a) (or any other 

suitable definition) and E^22 composite transverse modulus. For the linear 

case, the following relation holds 


®Z22T “ ^Z22p®Z22 


(7) 
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(a) Matrix. (b) Boron -epoxy composite. 

Figure Z -Matrix and composite stress strain curves. 


where 22p measured at the first knee (or point of linear deviation on the composite 
stress -strain curve as is illustrated in fig. 2(b)) . (These data were obtained under Air 
Force contract AFML-TR-66-313.) It is shown in chapter 2 of reference 1 that the 


strains and 


^I22p related by the following equation 


•I22p ~ P22T 


mpT 

li22 


( 8 ) 


where ^22T is the theory -experiment correlation factor and and <P ^22 ^.re the void 
effect and the matrix-strain-magnification factor, respectively, and are defined in the 
appendix. Substitution of equation (8) in equation (7) yields 


®Z22T “ ^22T 22 

li22 

Equation (9) relates S^22T limiting-matrix tensile strain, to the void effects, to 

the composite transverse modulus, and to the fabrication process through ^22T' 
interesting to note that both local effects (/3^ and <p and average effects (E^22^ 
fluence the transyerse tensile strength. The coefficient i322'p should be selected so 
that equation (9) correlates with experimental data 
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Transverse Compressive Failure Stress Sj22c 


The governing equation for the transverse compressive strength is derived in a 
fashion similar to that used for ^i22T!‘ result is 


c _ o ^mpC 

^122C ~ P22C — 

^v^li22 


( 10 ) 


where ^22C theory -experiment correlation factor, is the limiting-matrix 

compressive strain, and the remaining variables are the same as for 22T’ should 
be noted that equations (9) and (10) differ only in the correlation coefficients and the 
allowable matrix strains. 


Intralaminar Shear Failure Stress S 


I12S 


The governing equation for the intralaminar shear strength is derived by a procedure 
similar to that used for S^22T ‘ result is 


^Z12S ^12S 


mpS 
^v*^ jul2 


'U2 


( 11 ) 


where jSj 2 g is the theory -experiment correlation factor, is the allowable matrix 

shear strain, i® composite shear modulus, and /3^ and <P^j^2 void 

ans shear matrix-strain-magnification factors and are defined in the appendix. The void 
effects are the same for S^22T’ ®f22C’ ®Z12S (^®^- ^) • 


Selection of the Correlation Coefficients for Simple Strengths 


The theory-experiment correlation coefficients in equations (2), (3), (9), (10), and 
(11) are selected from simple experimental setups as follows (see fig. 1): 
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Coefficient 

Stress 


^122 

°"Z12 

%T’ ^mT 

®^11T 

0 


3 

^fC’ ^mC 

®Z11C 

0 



^22T 

0 

®Z22T 



^22C 

0 

^122C 



^12S 

0 

0 

®Z12S 


It is recommended that /3j2S evaluated from thin tubular test specimens rather than 
short beam specimens since the uniform shear strength of a ply is required. A numerical 
example in selecting is illustrated in the appendix. 

The important points to be noted in connection with equations (9) to (1 1) are that 

(1) The failure stresses S^22T’ ^Z22C’ ®Z12S sensitive to the matrix 

properties, to the composite elastic properties, to the void effects, and to the matrix - 
strain-magnification factor. Therefore, it is important that the matrix-strain- 
magnification factors be determined with sufficient accuracy. 

(2) The nearness of the coefficients ^22T’ ^22C’ ^12S unity is a measure 

of the validity of the hypothesis and of the insensitivity of the failure stresses S^22T’ 
S^22C’ ®Z12S to th® lubrication process. 

It should be clear from the discussion to this point that all five UFC failure stresses, 
filament and matrix properties, filament and void content, matrix- strain -magnification 
factors and the UFC elastic properties are needed to evaluate the correlation coefficients 
Thus, the failure stress as well as filament and matrix properties and filament and void 
content should be made available by the material supplier. It cannot be overemphasized 
that these properties need to be known accurately for meaningful formulations of failure 
criteria and in particular for the selection of theory-experiment correlation factors. The 
matrix-strain-magnification factors and the UFC elastic properties can be computed when 
the constitutent material properties are known (appendix and refs. 1, 23, and 24). 

The simple strengths of UFC from several filament matrix systems are available in 
the literature and are listed in table I. Three points should be noted in this table: (1) the 
simple strengths are for one filament volume content (k^; (2) the simple strengths of the 
last four composites are preliminary data and may be modified as more published data 
become available; and (3) the magnitudes of the tensile and compressive strengths are 
considerably different. The correlation coefficients selected from the UFC simple 
strengths in table I are listed in table n. In table II many of the correlation coefficients 
are near unity. Some exceptions are the coefficients for all the composites and 

^22T some Thornel composites. The point to be noted is that the rule of mixtures 
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TABLE I. - UNIDIRECTIONAL FILAMENTARY COMPOSITE 


SIMPLE STRENGTHS 


Material 

Failure or limit stress, ksi 

Filament 

content, 

kf 

Refer- 

ence 

®Z11T 

®aic 

®f22T 

®Z22C 

^fl2S 

Boron 

^210 

195 

8.1 

26.4 

12.1 

0.50 

20 

E -Glass 

157 

101 

*^4.0 

20.0 

•^6.0 

.50 

15 

S -Glass 

268 

207 

3.3 

21.0 

5.5 

.60 

30 

Thornel-25 

92 

67 

1.0 

21.0 

4.0 

.50 

27 

Thornel-40 

140 

91 

1.0 

19.0 

3.7 

.67 

27 

Thornel-50^’^ 

115 

60 

3.6 

17.0 

2.6 

.60 

31 

Morganite-I 

130 

120 

6.0 

o 

d 

8.0 

.50 

32 

treated^ ’ ^ 








Morganite-n 

150 

130 

8.0 

20.0 

11.0 

.50 

32 

treated^’ ^ 








Beryllium^ 

1 67.4 

i 60.0 

5.0 

22.0 

11.7 

: .50 

29 










^Different k^; Boron, 0. 54. . 

'^Estimated. 

^Estimated from data reported in the reference. 

*^These values are preliminary and might change as more information 
becomes available. 


(matrix strength limited) is not a representative mechanistic model for The re- 

sults of the Thornel composites indicate poor constituent bond, which behavior is similar 
to that of a large void content composite. The important point is that the equations de- 
rived here describe the UFC simple strength behavior satisfactorily. An additional im- 
portant point concerning the first-level semiempirical theory is that the correlations 
coefficients in table n can be used to analyze, design, and develop failure envelopes for 
UFC of filament and void contents in the practical range 0. 35 ^ kj ^ 0. 75 and 0 ^ k.^ < 

0. 20. This should hold so long as the fabrication process variables noted in equation (1) 
remain invariant for that particular process. Of course, if any of these fabrication pro- 
cess variables change, then the coefficients need to be reevaluated. 

The coefficients in the last four lines of table n are needed in the combined -stress 
strength criterion described in the next section. 
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TABLE II. - CORRELATION COEFFICIENTS FOR VARIOUS COMPOSITES 


[ The )3 coefficients correspond to data in table I. ] 


Correlation 

coefficient 

Equa- 

tion 

Fiber -matrix composite 

E -Glass 

S-Glass 

Boron 

Thornel-25 

Thornel-40 

Thornel-50 

Morganite-I 

Morganite-n 

Beryllium 

^£T 

(2) 

0 

.82 

1.00 

0.84 

1.00 

0.84 

0.83 


1.00 

0.84 

1.00 

^fC 

(3) 


.33 


.55 


.12 


.12 


.08 


.05 


.09 


.16 


.08 

^mT 

(2) 

1 

.00 

1.00 

1.00 

1.00 


1.00 

1.00 


1.00 

1.00 


L.OO 

^mC 

(3) 

1 

.00 

1.00 

.1.00 


L.OO 


1.00 

1.00 


1.00 

1.00 


l.OO 

^22T 

(9) 


.55 


.66 


,90 


.06 


.08 


.26 


.50 


.70 


.53 

^^220 

(10) 

1.10 

1.70 


1.22 


.49 


.50 


.50 


.65 


.70 


.90 

^12S 

(11) 


.86 


1.30 


1.50 


.48 


.46 


.27 


.91 

1.37 


1.40 


(4) 

13.30 

13.30 

13.3 

13.3 

13.3 

13'. 3 

13,3 

13.3 

13.3 


(4) 

31'900 

31 900 

31 900 

31 900 

31 900 

31 900 

31 900 

31 900 

31 900 

^Z12TT 

(12) 

1.00 

1.00 


1.00 


1.00 


1.00 


1.00 


1.00 

1.00 


1.00 

b 

^Z12CT 





















b 

^U2TC 



1 


















ir' ^ 

^H2CC 



) 


1 




1 













^For glass-epoxy composite; might need reevaluation for other composites. 
^These values were assumed; might need modification. 


COMBINED-STRESS STRENGTH CRITERION 


The governing equation for this criterion is derived from the following two postulates 
(1) at the onset of failure, the distortion energy under simple and combined loading re- 
mains invariant, and (2) the tensile and compressive properties of UFC are the same up 
to the onset of failure. These two postulates are based on the von Mises criterion for 
isotropic materials and on the experimental observation that the distortion energy of UFC 
remains invariant under rotational transformation (ref. 1, section 3.3). The formal 
derivations are described in detail in section 3. 3 of reference 1. The resulting equation 
(from ref. 1) is 

F(ct;,S;,K,,2).1- 

( 12 ) 


. k: 




722^7 


712Sy 




^n\a°l22^ 


1^122^ 


where F denotes the combined-stress strength criterion as follows: 


F(ct^,S^,K^12) > 0 

no failure 

F(a^,S^,K^12)=0 

onset of failure 

F(cr^,S^,K^j2) < 0 

failure condition exceeded 

and where subscripts a and /3 denote 
applied stress state determined from the 

T (tension) or C (compression), denotes the 

stress analysis, denotes the UFC simple 


strength either determined from the equations described previously or measured exper- 
imentally, and K^i2q!/ 3 theory-experiment correlation coefficient and is deter- 

mined as will be described subsequently. The coefficient K^j 2 given by (from ref. 1, 


ch. 3) 


K 


Z12 


^122 ^ 12 ^ 


[®^ ll®^22 + ^112 ^ 113 ^ ^121 + ^ 12 ^'\ 


(13) 


where and denote UFC modulus and Poisson's ratio, respectively. The sub- 
scripts 11, 13, etc. , refer to the corresponding axes in figure 1. For the case of iso- 
tropic material, equation (13) reduces to unity, as can be verified by direct substitution, 
and equation (12) reduces to the well known von Mises criterion. 
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Several important points should be noted at this juncture: (1) The UFC do not exhibit 
similar properties in tension and compression, as was stated in the second postulate 
made at the beginning of this section. Therefore, K^]^2 q!/ 3 introduced to compensate 
for this disparity. (2) Equation (12) describes failure at each quadrant by using, at most, 
four parameters. (3) The correlation coefficient K^j2a!/3 different values in 

different quadrants. (4) The product defined as one constant and 

determined experimentally. However, this disguises the composite effects which are 
introduced into equation (12) through (5) The product restricted 

to any range, that is, -°°"^ii2a^^ll2 Equation (12) is applicable 

to all materials exhibiting generally orthotropic elastic symmetry and is not restricted 
only to UFC. 

Values of the variable K^j2 against are plotted in figure 3 for several filament- 
resin systems. The graphs in figure 3 are applicable to composites with various void 



Figures. -Combined-stress strength -criterion 
coefficient (eq. (13)). 


contents since K^j 2 i® slightly sensitive to the void content. A procedure to select 
the coefficients K^i2aj3 illustrated in figure 4. The experimental results in this 
figure are for JT-50 graphite at 10 percent porosity as reported in reference 21. For 
this material K^j2 “ 0-85. As can be seen from figure 4, the coefficients 
be selected so that a good theory -experiment correlation can be obtained. One important 
point to be noted in figure 4 is that the slope of F(a^, S^, K^j2^ discontinuous across the 
quadrant junctures. This type of behavior is typical for orthotropic materials. Another 
important point is that the failure criterion is sensitive to K^i2q'/ 3 tension-tension 

(TT) and compression-compression (CC) quadrants. 
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o Experimental (ref. 21) 



Figure 4. - Evaluation of correlation coefficient K(i 2 op (eq. {12)). 


APPLICATIONS, RESULTS, AND DISCUSSION 

The two-level semiempirical failure theory can be used in several ways: composite 
failure analysis, design, structural synthesis, generation of strength envelopes, and as 
an aid to experimental work. Here, the discussion is restricted to composite failure 
analysis (which serves as a verification of the two -level theory) and to the generation of 
strength envelopes (which points out problem areas in testing and design). 


Composite Failure Analysis and Theory Verification 

The two-level, semiempirical failure theory has been applied to several multilayered 
filament-matrix composites. The input data in these applications consisted of the con- 
stituent material properties, the filament and void contents, the correlation coefficients 
(table n), the composite geometry and the failure or maximum load. The generation of 
other required properties (ply elastic constants, simple strengths, etc.), the composite 
stress analysis, and the failure test according to equation (12) were carried out by a 
multilayered-filamentary-composite-analysis computer code (ref. 25). Typical results 
are presented in tables in to VI and described subsequently. 

The results in table HI are for Thornel-25 epoxy composite (for the narrow flat 
specimens reported in ref. 26) . The first four columns in table HI contain the composite 
geometry, the fifth column the applied load, and the last two columns the value of the 
failure criterion (eq. (12)) for the first and second plies. As can be seen in the last two 
columns of table HI, the criterion predicts failure or nearly so in all composites except 
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TABLE m. -RESULTS OF FAILURE ANALYSES OF THORNEL- 


25-EPOXY MULTILAYERED COMPOSITES BASED ON 


NARROW SPECIMEN DATA FROM REFERENCE 26 


Composite 

Fiber 

Ply 

Failure 

Failure criterion 

ply 

content, 

thickness, 

load, 

(eq. (12)), 

arrangement 




F(<r^,S 

1 ’ 12^ 



in. 

Ibf/in. 







First ply 

Second ply 

9(0) 

0.324 

0.0158 

9270 

0.16 

0.16 

9(0,90) 

.524 

.0136 

2990 

.85 

-1.25 

9(90, 0) 

.557 

.0131 

3650 

-2.87 

.52 

9(±45) 

.592 

.0127 

876 

-.24 

-.24 

9(±45) 

.607 

.0126 

822 

-.12 

-. 12 

4(0,90,90,0) 

.507 

.0138 

2170 

.49 

^-41.6 

4(90,0,0,90) 

.475 

.0142 

1922 

^-32.2 

.60 

3(0,90,0) 

.538 

.0134 

1726 

.69 

^-28.7 

3(90,0,90) 

.554 

.0131 

753 

^-16.9 

.57 

3(±45) 

.428 

.0150 

^189 

.45 

.61 


^Load transfer difficulties. 
^No failure 


TABLE IV. - RESULTS OF FAILURE ANALYSES OF THORNEL-40-EPOXY 


MULTILAYERED COMPOSITES BASED ON DATA FROM REFERENCE 27 


Composite ply 

Fiber 

Ply 

Failure load, 

Failure criterion 

arrangement 

content, 

thickness, 

Ibf/in. 

(eq. 

(12)). 


^f 

h' 

N 

N 




in. 

X 

y 

First ply or 

Second ply or 






first pair 

second pair 

4(90,0.0,90) 

0.64 

0.0086 

1960 

0 

®^-16.5 

0.63 




0 

2300 

.49 

^-23. 1 




-1520 

0 

.97 

-.03 




0 

-1750 

-.36 

.95 

4(90,10,-10,90) 

0. 57 

0.010 

1440 

0 

^-7.80 

0.73 




0 

2240 

.57 

^-17.0 




-1520 

0 

.96 

.14 




0 

-1960 

.25 

.88 

4(10,-10,-10,10) 

0.69 

0.0079 

2120 

0. 

^0.52 

0. 52 




0 

31.6 

.23 

.23 




-1710 

0 

-2.18 

-2.18 




0 

-15. 5 

-.16 

-.16 

6(10, -10,45. -45, 10. -10) 

0.56 

0.0102 

3300 

0 

0.30 

-0.67 




0 

428 

-4.80 

-.78 




-2940 

0 

.65 

-1.09 


^Load transfer difficulties 



TABLE V. - RESULTS OF FAILURE ANALYSES OF E -GLASS-EPOXY COMPOSITES 


BASED ON LONG TUBULAR SPECIMEN DATA FROM REFERENCE 28 


Cylinder^ 

Composite ply 
arrangement 

Fiber 

content, 

kf 

Ply 

thickness, 

h’ 

in. 

Failure 

load, 

Nx> 

Ibf/in. 

Failure criterion 
(eq. (12)), 

+9 or -9 

90° 

4 

12(45, -45,90,90,...) 

0.632 

0.00455 

-568 

-0.25 

0.80 

5 


.625 

.00461 

-597 

-.30 

.79 

6 


.655 

.00372 

-496 

-.57 

.75 

7 

12(67.5,-67.5,90,90,...) 

0.663 

0.00464 

-541 

1 

o 

o 

CD 

CO 

CD 

d 

8 


.660 

.00450 

-491 

.06 

.72 

9 


.694 

. 00460 

-557 

-.32 

.61 

1 

12(25,-25,90,90,...) 

0.638 

0.00459 

-700 

0.19 

0.83 

2 


.650 

. 00457 

-654 ■ 

.22 

.84 


lumbers refer to cylinders reported in table III of reference 28. 


TABLE VI. - RESULTS OF FAILURE ANALYSES OF BERYLUUM 
MULTILAYERED COMPOSITES BASED ON DATA 


FROM REFERENCE 29 


Composite ply 

Fiber 

Ply 

Failure 

Failure criterion 

arrangement 

content, 

thickness, 

load, 

(eq. 

(12)), 


^f 


N , 

X’ 

Ibf/in. 




in. 

First ply 

Second ply 

6(0,90,0,0,90,0) 

0.443 

0.00666 

2693 

-0.04 

0.06 


.475 

. 00643 

2844 

-.16 

.02 


.541 

. 00602 

3012 

-.10 

-.53 

6(0,90,0,0,90,0) 

0.443 

0.00666 

-2565 

“0.23 

0.99 


.475 

. 00643 

-2732 

^-.14 

.99 


.541 

.00602 

-2997 

®.31 

.99 

6(0,60,-60,-60,60,0) 

0.459 

0.00654 

1787 

0.12 

0.57 


.459 

. 00654 

-1915 

^.38 

.95 


^Based on estimated value for 




the last one which was not carried to failure . The large negative values of the criterion 
indicate primarily two possibilities: (1) Load was nonuniformly transferred from the 
outer to the inner plies . (2) The transverse plies failed at an early state of the loading 
process and the load was carried primarily by the longitudinal plies. No stress-strain 
graphs were recorded to identify this condition. The criterion values of the longitudinal 
plies are in accord with the second possibility. This is apother important use of the 
semiempirical theory; that is, it points out problem areas which need be either remedied 
or avoided. The results presented in table IV are for Thornel -40 -epoxy composites (flat 
specimens) . The experimental data are reported in reference 27 . The criterion predicts 
failure (or nearly so) for all these composites. The explanation of the large negative 
criterion values in this table is the same as that for table m. 

Results of a different type are presented in table V. The experimental results for 
these cases were used to evaluate buckling of E -glass-epoxy cylinders and are reported in 
table in of reference 28 (Cylinders 1, 2 and 4 to 9). The nmnber of plies and the ply 
thicknesses were deduced from the data presented in this reference. Here, the combined- 
stress strength criterion is applied to determine possible ply failure at the reported buck- 
ling load. The results of the analysis indicate failure of the ply material for the ±45° 
plies in cylinders 4 to 6 (reported to have failed by buckling) and the ±67. 5° plies in cy- 
linders 7 to 9 (reported to have failed by buckling and material failure) and no material 
failure, but nearly so, for the ±25° plies in cylinders 1 and 2 (reported to have failed by 
buckling). The results of the analysis of these cases further illustrate the usefulness of 
the criterion proposed herein in interpreting experimental results and also the failure 
mode complexities of angle ply composites. 

Results of analyses for beryllium -epoxy composites (flat specimens) are presented 
in table VI. The experimental data for these composites were reported in reference 29 . 
The ply -thickness simple strengths and the failure load (failure stress only reported) 
were deduced from the data presented in reference 29 . The criterion predicts failure 
for the tensile load cases but for only one of the compressive load cases. This is be- 
cause of the difficulty encountered in establishing a reasonable longitudinal compressive 
failure stress from the data reported in reference 29. 

The results presented in tables HI to VT tend to substantiate the two -level theory 
proposed here. And what is more important, they illustrate a relatively simple pro- 
cedure to predict composite strength behaviors from basic constituent properties and 
fabrication process considerations. 

Generation of Strength Envelopes 

Envelopes for simple strengths and for the combined-stress strength criterion are 
presented in this section. Envelopes for unidirectional composites loaded with normal 
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Fiber content, kf cs-49311 


Figure 5 . - Limit stress for S-glass-epoxy compos- 
ites (based on correlation coefficients in table II 
and zero void content). 



Figure 6 . - Limit stresses for Thornel -50 - epoxy com- 
posites (based on correlation coefficients in table II 
and zero void content). 
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Figure 7. - Limit stresses for boron -epoxy compos- 
ites (based on correlation coefficients in table II 
and zero void content). 


and shear loads at some angle from the filament direction are also presented. 

Figures 5 to 7 present graphs of the simple strengths 
(eq. (4)), S^22T’ ®^22C’ ®U2S filament content kj with zero voids. These 

figures are, respectively, for S-glass, Thornel-50, and boron -filament UFC. The 
graphs in these figures were generated from the correlation coefficients in table I. The 
important points to be noted from these figures are that (1) the transverse and shear 
limit stresses decrease with increasing filament content for S-glass and boron compos- 
ites but remain rather invariant for Thornel-50 composites, (2) the decrease of these 
limit stresses S^22T’ ^122C’ rapid for boron composites at high 

filament volume content values, and (3) test results for transverse and shear properties 
for isotropic filament composites should be reported with accurate volume content par- 
ticularly in the high range. Those for orthotropic filaments (Efn/E£22 
be very accurate. If longitudinal compressive failure of UFC is governed by j jp 



Figure 8. - Limit stresses for Thornel-50 - epoxy 
composites (based on correlation coefficients in 
table II and 5 percent void content). 


20 



(constituent debonding and intralaminar shear), then this failure strength decreases with 
filament content (The value of kj should be measured fairly accurately) . 

Figures 8 and 9 present the simple strength envelopes for Thornel-50 and boron 
composites with 5-percent void content, respectively. Superposition of figure 6 with 8 
and figure 7 with 9 reveals a considerable drop in the transverse and shear strengths of 
composites with voids. It is important, therefore, to report the void content accurately 
when presenting experimental results on transverse and shear strengths. Of course die 



Figure 9. - Limit stresses for boron-epoxy compos- 
ites (based on correlation coefficients in table II 
and 5 percent void content). 


void content should be reported with longitudinal compressive strength as well since this 
strength could be governed by (constituent debonding and intralaminar shear). 

The combined-stress strength behavior for a Thornel-50 epoxy UFC is illustrated in 
figure 10. The contours in this figure represent strength envelopes for various values of 




to 


CS-49312 



Figure 10. - Combined -stress strength -criterion for Thornel-50 - 
epoxy composite (based on data in table I and eq. (12)). 
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intralaminary shear expressed as a fraction of three important points to be 

noted in figure 10 are that (1) Thornel-50 composites under the proper proportion of 
combined loading can resist considerably more load than their simple strengths would 
indicate, (2) the normal load capacity of this UFC is insensitive to small values of shear 
loads, and (3) the longitudinal compressive strength of UFC is very sensitive to 

transverse tensile loads while the longitudinal tensile strength is sensitive to transverse 
compressive loads and very sensitive to small transverse tensile stresses. Therefore, 
it is important to bear in mind (when testing for longitudinal strengths) that the trans- 
verse stresses should be completely eliminated. A small amount of shear stress can be 
tolerated. 

Figure 11 illustrates the strength envelope of UFC when loaded with normal loads at 
some angle to the filament direction. The upper part of the figure is for tensile load and 
the lower for compressive. These envelopes were obtained by expressing cr^ in equa- 



Figure 11. - Limit stress for off-axes normal 
load boron -epoxy composite. (Limit stresses 
from fig. 7 at 50 percent fiber content. ) 


tion (12) in terms of and then solving the resulting expression for a . As can be 
seen in this figure, the strength drops off very rapidly. The composite has lost about 
50 percent of its strength when the load is applied 5° to the filament direction. It is 
imperative, therefore, to have the load completely alined with the filament direction 
for longitudinal tensile strength tests. In tests to determine strengths for 6 > 40° (also 
transverse strength 6 = 90°) the load alinement is not critical. 

Figure 12 Illustrates an analogous effect for the case of shear load. Both positive 
(tending to elongate the filaments) and negative shear curves are shown in the figure. 

One important effect brought out by the negative shear curve (not widely recognized) is 
that UFC loaded with negative shear load 45° to the filament direction has a very low 
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Figure 12. - Limit stress for off-axes shear load 
boron-epoxy composite. (Limit stresses from 
fig, 7 at 50 percent fiber content. ) 


shear strength. This means that, if ±45® plies are introduced to help carry shear loads, 
the -45® ply fails (perhaps not completely) at relatively small loads and the ±45° ply 
carries the load. The partially failed -45® ply causes the composite to exhibit a non- 
linear load response as the loading increases. Other important points to be kept in mind 
for shear strength tests and design are that (1) the shear strength of UFC is insensitive 
to small angular deviations (0 0 and 0 90 in fig. 12), (2) the strength of UFC loaded 

with positive shear load at 45® is sensitive to small angular deviations, but it is insen- 
sitive if loaded with negative shear load, and (3) the strength of a UFC loaded with pos- 
itive shear load at 45® is approximately three times greater than the similar case with 
negative shear load. This ratio is approximately equal to 22C^®Z 22T’ 

Several recommendations have already been made on how theoretical work can aid 
the experimental effort. On the other hand, controlled experimental work is a very 
important asset in formulating and verifying any theory. For this reason, it is impera- 
tive that the complete test record (type of specimen, constituents, voids, type of test 
and techniques, type of failure, strain rate, means for measuring strain and elongation, 
and any other factors which influence the result) be reported when presenting experimen- 
tal data for filamentary composite properties. 


CONCLUSIONS 

A two-level, semiempirical theory was developed to predict the strength behavior of 
unidirectional filamentary composites (UFC) from constituent material properties and 
from fabrication process considerations. This two-level theory describes the strength 
behavior of several filament matrix composites reasonably well. The two-level theory 
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can be used to generate simple and combined strength envelopes for UFC. 

The results of this investigation lead to the following conclusions: 

1. The simple strengths are sensitive to the correlation coefficients and thereby 
to the particular fabrication process. 

2. The simple strengths are fairly sensitive to void content. 

3. The transverse and shear strengths are decreasing functions of the filament con- 
tent (particularly at the high range) for isotropic filament composites. 

4. The combined -stress strength criterion is sensitive to its correlation coeffi- 
cient in the tension-tension and compression-compression quadrants and relatively 
insensitive in the other two quadrants. 

5. The longitudinal compressive strength is very sensitive to the presence of trans- 
verse tensile stress. 

6. The longitudinal tensile strength is very sensitive to small transverse compres- 
sive stress. 

7. The normal load carrying capacity of a unidirectional filamentary composite 
decreases rapidly as the angle between filament and load direction increases. Con- 
sequently, the longitudinal tensile and compressive strengths are very sensitive to load 
misalinement. 

8. The shear strength is insensitive to load misalinement. 

9. The positive shear load carrying capacity of a unidirectional filamentary com- 
posite increases as the angle between filament and load direction increases to a max- 
imum at 45°, but it decreases for negative shear load to a minimum at 45°. 

10. Plies introduced to resist shear forces may do so by nonlinear response. 

11. The complete test record should be reported when presenting experimental data 
on composite materials. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 23, 1969, 

124-08-06-01-22. 
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APPENDIX - USEFUL RELATIONS 


Actual Filament and Matrix Volume Content 

Let kj and denote the apparent filament and matrix volume content, respec- 
tively, and let k^, kf, and k^^ denote the actual void, filament, and matrix volume con- 
tent, respectively. Then it can be shown (ref. 1, appendix A) that 


m 


(1 - 

(Al) 

- k^)(l - kf) 

(A2) 

= 1 -kf 

(A3) 


Matrix-Strain-Magnificatlon Factors and Void Effects 

The transverse and shear matrix-strain-magnification factors are, respectively, 
given by (ref. 1, ch. 2): 


(p 


p.22 


1 + p(A = 1) 


and if E^j^o-^22 ^a2^Z22^Zll " ° 


1+P(^fl2 - ^12^1 


'^ 122^111 " ^Z 21^ai^^22 
®ai‘^Z22 " '^Z12^Z22^Z11 


‘^M22 


= 1 


(A4) 


(A5) 


*^^12 


1 - p 1 


^ml2 l 

%2, 


(A6) 


The maximum void effect factor is given by (ref. 1, ch. 2): 
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where 


P = 



(A8) 


__ 1 - ^m22 

^ “ % 12%21 %22 


(A9) 


and E, G, v, and a denote longitudinal modulus, shear modulus, Poisson's ratio, and 
stress, respectively. The subscripts v, f, and m denote void, filament, and matrix, 
respectively, and the numerical subscripts correspond to the filament directions depicted 
in figure 1. The variable is the void content, and kj and k^ are defined by equa- 
tions (Al) and (A3), respectively. 


Evaluation of Correlation Coefficients for Simple Strengths 

Several ways can be used to evaluate the correlation coefficients. The simplest one 
is illustrated by the following example using Thornel-25 epoxy: 

= 180 000 

IIT “ 
kj = 0. 50 

k^ = 0. 5 
m 

Efii = 25x10® psi 
Emil = 0. 55x10® psi 



Solving equation (2) for yields 





^mll\ 

%11 / 


(AlO) 


Assuming = 1.0 and substituting the values for all the variables in equation (AlO) 
yield <^1.0. The remaining correlation coefficients can be evaluated in the same 
fashion. 
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